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ABSTRACT 

We consider the existence of Hamiltonians for autonomous non-holonomic mechanical systems. 
The approach is elementary in that the existence of a Hamiltonian for a non-holonomic system 
is equivalent to the existence of an appropriate Lagrangian for the system in question. The 
existence of such a Lagrangian is related to the inverse problem of constructing a Lagrangian 
from the equations of motion. A simple example in three dimensions with one non-holonomic 
constraint is analyzed in detail. In this case there is no Lagrangian reproducing the equations 
of motion in three dimensions. Thus the system does not admit a variational formulation in 
three dimensions. However, the system in question is equivalent to a two-dimensional system 
which does admit a variational formulation. Two distinct Lagrangians and their corresponding 
Hamiltonians are constructed explicitly for this two-dimensional system. 
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1 Introduction 



Hamilton's principle for mechanical systems with non-holonomic constraints has recently been 
discussed by Flannery pQ. In particular a variational formulation of the equations of motion 
of a mechanical system was discussed both for holonomic and non-holonomic constraints. It 
was shown that while the equations of motion for a system with holonomic constraints can 
be obtained as variational equations, with the constraints being taken into account by the 
multiplication rule in the calculus of variations [2], the corresponding procedure with non- 
holonomic constraints leads to equations which differ from the correct equations of motion. 

The problems discussed by Flannery are not new; they have been discussed in the literature 
at least since Hertz's textbook |3|, in which the use of variational principles in mechanics was 
questioned. Two papers published by Jeffreys [J] and Pars [5] consider Hamilton's principle for 
non-holonomic systems, and propose rectification of previous papers in which the variational 
procedures discussed by Flannery had been used also for non-holonomic systems. 

Several papers have advocated the use of a variational principle involving the multiplication 
rule in the calculus of variations for non-holonomic systems. In addition to the papers of this 
kind quoted by Flannery and by Pars and Jeffreys, respectively, we mention a paper by Berezin 
[B], in which no distinction is made between holonomic and non-holonomic systems. 

It is also appropriate to to mention that, in contradistinction to the original "Classical 
Mechanics" by Goldstein [7] , the 3rd edition of this classical mechanics textbook advocates the 
use of a variational principle involving the the multiplication rule for non-holonomic systems 
[8]. However, the use of this principle for non-holonomic systems was later retracted [9]. This 
fact was pointed out already by Flannery. 

It appears that if a system with non-holonomic constraints does not admit a variational 
formulation, then the dynamics of the system is not governed by a Hamiltonian H . This is 
the question we address in this paper: Can a non-holonomic system be described in terms of 
Hamiltonian equations of motion? We confine the detailed discussion to a simple example in 
three dimensions introduced by Pars [5] . We show that in this case the equations of motion are 
reducible to a set of equations for a two-dimensional autonomous system, which can be formu- 
lated as Hamiltonian equations. However, the original equations of motion in three dimensions 
do not admit a Hamiltonian formulation. 

Our analysis is elementary in that the existence of a Hamiltonian for a given non-holonomic 
system is considered to be equivalent to the existence of an appropriate Lagrangian L(q, q) for 
the system in question. By appropriate is meant that the Lagrangian is non-degenerate, i.e. 
that the equations defining the canonical momenta p 7 , 



are solvable for the generalized velocities <jj. It should be noted that we discuss only autonomous 
systems. Hence the Lagrangian is allowed to depend on time only through the coordinates q 
and velocities q. 
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The existence of an appropriate Lagrangian is related to the inverse problem of constructing 
a Lagrangian from the appropriate equations of motion. To the best of our knowledge, a 
complete solution to the inverse problem is not known in the general n-dimensional case for 
n > 3. 

In the next section we consider the Lagrange equations of motion for an autonomous me- 
chanical system with both holonomic and non-holonomic constraints. This problem was con- 
sidered in some detail by Flannery [lj. For the sake of completeness we consider the equations 
obtained from the generalized form of d'Alembert's principle and the equations which follow 
from a variational procedure with constraints implemented by the multiplication rule. In the 
non-holonomic case these equations are not identical. 



2 Lagrange equations with constraints 

Consider an autonomous mechanical system with independent generalized coordinates qi, q n , 
and velocities tji, q n . Let the kinetic energy be T, and the generalized applied forces on the 
system be Q^,j = 1, ...,n. The generalized principle of d'Alembert (see e.g. the classical texts 
by Goldstein [7J or Whittaker [10]) then gives the following equation, 

£{i®-S-*}*-* 

where the quantities Sqj are virtual displacements of the system. If the virtual displacements 
Sqj,j = 1, . . . ,n are independent, then the equation fl2]) results in the Lagrange equations of 
motion, 

d_ fdT\ _dT =QJj = 1 n (3) 
dt \dq~j J dqj 

We generalize to systems with 1 < m < n independent non-holonomic constraints, which 
are taken to be linear and homogeneous in the velocities. The constraint equations are of the 
following form, 

n 

J2 a i 3 (l^ ->?n)?7 = 0, % = 1, ...,m < n, (4) 

3=1 

where the quantities a/ ,i — 1, . . . ,m, j = 1, . . . ,n, are given functions of the variables qi, q n . 

The derivation given below of the equations of motion for this non-holonomic system can 
be found e.g. in the textbook by Whittaker [TO] . 

Implement the constraints i^j by regarding the system to be acted on by external applied 
forces Qi and by certain additional forces of constraint Q'\ j = 1, . . . , n, which force the system 
to satisfy the non-holonomic conditions (Tj0). The equation (EJ) is then replaced by the following 
equation, 
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In Eq. (J5} the virtual displacements Sqj, j = l,...,n, can now be regarded as independent. 
Thus one obtains the equations of motion, 

d (dT\ dT , . 1 

The forces of constraint, Q'\ j = 1, . . . ,n, are a priori unknown, but they are such that, in 
any instantaneous displacement 5q> ,j = 1, . . . ,n, consistent with the constraints (j3|, they do 
no work. The non-holonomic constraints (j3J) imply the following conditions on the possible 
instantaneous displacements 5q°,j — 1, . . . ,n of the system, 

n 

J2ai J (qi,-,qn)Sqj = 0, i = l,...,m < n. (7) 
i=i 

For any instantaneous displacements Sq^j = l,...,n, which satisfy the conditions ([7]), the 
work 5W done by the constraint forces Q'- 7 , j = 1, . . . ,n equals zero, 

n 

sw = <3'W = o. (8) 

3=1 

The conditions (0) and (jHJ) together imply that 

Qtf = ]TAV,i = i, ( 9 ) 

i=l 

where the quantities X\ i = 1, . . . , m, are time-dependent parameters. The equations ([6]) have 
been reduced to 

d ( 8T\ dT m 

s©-fe =<J,+ £ AV>i=1 '" >B - (10) 

To these n equations of motion one should add the m equations of constraint (j3J). We have n+m 
equations for the determination of n + m quantities qj(t),j — 1, . . . ,n, and A l (t), i = 1, . . . , m. 

It should be observed that in the argument above, one has not required the constraint 
equations (jlj) to be in force under general variations qj — > qj + Sq^; the constraints (j4]) are only 
imposed on the actual motion of the system. 

Now assume that the external applied forces Q 3 ,j = 1, . . . , n, can be expressed in terms of 
a potential V such that, 

Q J = ~— + -(— ) 3 = 1 n (11) 
dqj dt \d(jj J 

Using the notation 

L :=T-V, (12) 

the equations (flOj) can be written as, 

d I dL \ 3L ™ ( t . 



It should be observed that the m one-forms ([7]) are non-integrable by assumption, for oth- 
erwise the system would be holonomic. In the integrable case (after multiplying the conditions 
(j4j) with integrating factors if necessary) one would have 

a/ = -—,« = l,...,m, (14) 
% 

where the functions G{,,i = 1, ...,m, are m independent functions of the variables qj,j = 
l,...,n, 

G l = G i (q 1 ,...,q n ),i = l,...,m. (15) 

The m constraint equations (jlj) would then be equivalent to the following m holonomic con- 
straints, 

Gi(qx, . . . ,q n ) = d,i = 1, . . . ,m, (16) 

where the quantities Ci,i = l,...,m are constants. In this case the equations ( fl3l) are the 
Euler-Lagrange equations of the variational problem 

S J dtL = 0. (17) 

under the constraints fflBT) . These constraints can be implemented with the multiplication rule 
in the calculus of variations. This leads to the following free variational problem with Lagrange 
multipliers X\ i = 1, . . . , m, 



5 / dt 



Lq + ^ A* {Gi — Ci 



i=l 



0. 



The variational problem yields 

d (8L \ dL f,.,9G, . 

*l^j-^--§ A ^ = ' J=1 "'"' i ' <19) 

The system of equations (ITP1) . together with the constraints ( TTB|) . are the correct equations of 
motion for the system under consideration in the integrable (holonomic) case. These equations 
are a set of Euler-Lagrange equations with the integrand in Eq. (Tl8|) as a Lagrangian L, 

m 

L:=L + y £\ i (G i -C i ), (20) 

i=i 

provided one adjoins the time-dependent parameters A l ,i = l,...,m, as new coordinates to 
the system. It should be noted that the usual method of transition to a Hamiltonian from the 
Lagrangian ([20]) does not apply, since the momenta conjugate to the new coordinates A J vanish 
identically. 

Contrary to the assertions in some of the papers referred to in the references p], [I], [5], 
as well as in reference [6], a similar procedure in the non- holonomic case does not lead to the 
correct equations of motion. Specifically, if one considers the variational problem ffl7|) under 



4 



the constraints (j3J) using the multiplication rule, one is led to the following free variational 
problem, 



5 / dt 



m n 



i=l i=l 

where the Lagrange multipliers are now denoted by fi l , i = 1, . 
following from Eq. (j2ip are, 



m 



= 0, (21) 
, m. The variational equations 





da/ 
d qj 



qk 



0. 



(22) 



The equations (1221) are not identical to the correct equations of motion ([TBI for the non- 
holonomic system under consideration. However, if the integrability conditions 



da,; 



da, 



dq k dqj 



0,i 



m, j, k 



, n. 



(23) 



are valid, in which case the system becomes holonomic, the equations of motion (122"]) coincide 
with the corresponding correct equations of motion ( ITBl) [equivalently Eqns, ( 1T91) ] upon a change 
of notation fi l — > A 1 , i = 1, . . . , m. 

The generalized principle of d'Alembert differs from the variational principle involving the 
multiplication rule in the case of non-holonomic constraints. One consequence of this difference 
is the fact that the equations of motion following from the principle of d'Alembert differ in form 
from the equations of motion which follow from the variational principle. It is not excluded that 
these equations may have the same solutions, however. It thus remains to consider the question 
of whether the equations of motion (1T3|) and the variational equations (122]) can have coinciding 
solutions in general. In his discussion of this problem Pars [5] used a simple yet non-trivial 
example in three-dimensional configuration space to show that the equations (1131) and (122]) in 
that case can not have coincident general solutions. We will give a detailed discussion of Pars' 
example below, adding a few details related to the relevance of the initial values. For clarity, we 
also pay attention to the dimensions of the quantities in the example, by including appropriate 
dimensional constants. 



3 Pars' example 

The example considered by Pars [5] is the case of an otherwise free particle of mass m in three- 
dimensional Euclidean space with coordinates designated by (x,y,z), except that the motion 
of the particle is subjected to the following non-holonomic constraint, 

zx-£y = 0, (24) 

where £ is a constant with the dimension length. The dimensional parameters m and £, which 
were absent in the formulation given by Pars, are introduced here for clarity. 
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The Lagrangian Lq in this case is the following, 

L = j(x 2 +f + z 2 ). (25) 

The Lagrange equations of motion ([TBI reduce now to the following equations, 

mi = Xz , my = —X£ , m'z = 0, (26) 

where A is a parameter, which is to be determined by solving the system of equations f]2"6"]) and 
(1231). The initial conditions £1X6 cLS follows, 



(x,y,z)\± =o = 0, (27) 

and 

(x,y,z)\ t=0 = (u, Q,w), (28) 

where u and u> are parameters at our disposal, except for the conditions u ^ and to^O. If 
u = or u> = then the solutions to the equations (|26|) with the initial conditions (1271) and 
( T28|) are trivial and uninteresting. It should also be noted that the condition y{0) = in ([28"]) 
above is not a free choice, but a consequence of the constraint equation ([21]) and the initial 
values (127]) . 

We now consider an alternative form of the Lagrange equations ( 1261) and the constraints 
(124"]) . Differentiating the constraint equation ([23D one obtains 

zx + xz - £y = 0. (29) 

Eliminating the quantities x, and y from the equation (I29p above with the aid of the equations 
([26]) . one obtains the following expression for the quantity A: 

mxz . . 

X = -¥T7- (30) 

Inserting the expression (15U]) for the parameter A into the original equations (12B1) one obtains 
the following three equations, 

%± = ~ 2 » ( 31 ) 

^ + Z" 1 

' r: (32) 



y £ 2 + ' 
and 

2 = 0. (33) 
It should be noted that the mass m does not appear in the equations ( |3T]) - ( 133] ) . 



Before proceeding further, we demonstrate that the set of equations (13TJ) - (1331 above are 
essentially equivalent to the original equations (126]) and the constraints (12^]) . Multiplying the 
equation ( 13TJ) with z, and subtracting the equation ( 132]) from the result, one obtains 

z'x - ty = -xz <^> - ^2/) = ^ (zi - iy) = C, (34) 
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where C is a constant. Using finally the initial conditions (l2~Tj) and (1281) to evaluate this constant 
one obtains 

C = 0. (35) 



The equations (1311) - (1331) thus imply the constraints ( 1241 when one also uses the information 
encoded in the initial conditions ( 1271 and (1281) . The equations ( 13~T1) - ( 1331) are indeed of the form 
(126|) . where the parameter A is identified with the expression (130|) . We have now demonstrated 
that the equations (l3Tj) - (1331) together with the initial conditions (l2~Tj) and (128|) are equivalent to 
the original equations (1261) and the constraints (1241) . The wording "essentially equivalent" used 
above was meant to reflect the fact that one had to invoke the initial conditions (12"T|) and (I28p 
in order to show that the constraints (1211) are a consequence of the alternative equations (13"T1) - 
and not a separate condition, as in the formulation (l2"B"j) which involves the parameter A. 

We note that the equation ( 13~T1) can be integrated, yielding 



xV£ 2 + z 2 = £u, (36) 

were u is the initial value at t = for the quantity x. The result ( 1361) will be used shortly. 
Consider then the variational problem (l2~Tj) for the case at hand, i.e. 

5 J dt [L - \i(zx - £y)\ = 0, (37) 



where the function L is given in Eq. (125|) . The differential equations which follow from Eq. 
(137|) are the following 

— [mx — \iz) = 0, (38) 
j t (my + fi£) = 0, (39) 

and 

m'z + fix = 0. (40) 
To the equations (138|) - (j40l) one should still add the constraint equation (pH|) . 

It will be shown that the solutions to the equations (1311) - (13"3"|) with the initial conditions 
(J27J) and (125|) can not satisfy the variational equations (13"%1) - (1411 and the constraint equation 
(124|) . except in certain trivial cases cases, as shown below. 



Assume now that there are appropriate solutions x(t), y(t) and z(t), which satisfy both sets 
of equations (T3Tj) - (1331) and (1381) - ( 1401) together with the constraint ( 1241 . respectively, under 
the initial conditions §7h and (ES]). From Eq. ([33]) and Eq. flU then follows that 



/ii = 0. (41) 
From the equations (jH] and (136|) then follows that 

//u = 0. (42) 
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The condition (jHJ, or equivalently (j42p . is thus necessary for the existence of functions x(t), 
y(t) and z(t), which satisfy both the set of equations (l3Tj) - (l33l) and the set of equations (138|) 
- f )40|) together with the constraint (l24"j) under the initial conditions fl27|) and (1281) . There are 
three possible cases to be considered: 

^ = 0, u = 0, (43) 
H = 0, w + 0, (44) 

and 

^ ^ 0, u = 0. (45) 

If the conditions (|43p are valid, then one one finds readily that the only common solutions of 
the equations (J5T]) - fl33|) and the equations (1381) - (T401) as well as the constraint equation (T24|) 
which satisfy the initial conditions are the following, 

a; = 0, y = 0, z = tyt. (46) 
Likewise, if the conditions (1441) are in force, then the only possible solutions are 

x = ut, y = 0, z = 0. (47) 
Finally, if the conditions ( 1451) are valid, one finds the following solution, 

x = 0, y = 0, z = 0, /i = c, (48) 

where c is a constant. 

The solutions ( 1461 ( HTj) . and ( 1481) . respectively, are the only functions which satisfy both the 
Lagrange equations of motion in the form (13T1) - (|33l) and the variational equations (1551) - ( j40l) 
together with the constraint (1241) . under the initial conditions (1271) and (128|) . These solutions 
are clearly exceptional in that the non-holonomic constraint (124"]) is no constraint at all for these 
solutions. 



4 Existence of Lagrangians and Hamiltonians in Pars' 
example 

It should be observed that the fact that the variational procedure involving the multiplication 
rule does not lead to equations of motion identical to those which follow from the generalized 
principle of d'Alembert in the case of non-holonomic systems, does not prove that there is no 
variational principle at all for non-holonomic systems. One may still wonder whether non- 
holonomic systems may nevertheless admit some kind of variational formulation. A straightfor- 
ward answer to this question is obtained if one can show that the correct equations of motion 
(fTBl) together with the constraints (T4|) constitute a set of Lagrangian equations with some ap- 
propriate Lagrangian. This is an inverse problem, which is trivial in the case of one- dimensional 
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systems. Complete results on the inverse problem in question exist for two-dimensional systems, 
but not for systems of dimension three or higher. We analyze the problem posed here only in 
the non-holonomic three-dimensional special case considered by Pars, which was analyzed in 
some detail above. 

The question is now whether the equations (13T1) - (1331) are the Euler-Lagrange equations with 
some appropriately chosen Lagrangian, or linearly equivalent to such Euler-Lagrange equations 
in three space dimensions. For this problem we refer to a paper by Douglas [TT] on the inverse 
problem in the calculus of variations as well as to a paper by Crampin et al. [12], which gives 
a geometric formulation of the inverse problem, with due reference to the paper of Douglas. 

Using results given in the papers by Douglas and Crampin et al, referred to above, one finds 
that the equations f[3"Tj) - (1551) can not be recast into linearly equivalent equations involving 
three variables, such that these equivalent equations are the Euler-Lagrange equations of some 
appropriate functional. We know that the space of dynamically accessible paths in the problem 
under consideration is in fact two-dimensional, so it is then natural to look for a variational 
formulation in a two-dimensional space. It will be shown that the system of equations (TBTj) 
- (J2SD can be reduced to an equivalent two-dimensional autonomous system, for which there 
exist Lagrangians. 



Eliminating the quantity x from equation in the system (1321) with the aid of the relation 
361) above, one obtains a two-dimensional autonomous system from the equations (13TT) - ( 1331) 
which involves the variables y and z only, 
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y = l 2 uz(f + z 2 ) % (49) 
z = 0. 

The simple system of equations (H9|) is indeed obtainable from a principle of stationary action 
in a space of two dimensions. There is in fact more than one Lagrangian for which the equations 
( 1491) are the Euler-Lagrange equations. It is known that Lagrangians which are derived from 
the equations of motion are not necessarily unique. [13] We display below two such Lagrangians 
Li and Ln, whose difference is not a time derivative of some appropriate function: 

Li := myz log — , (50) 



VF+z^ °Vc 

where Co is a constant with the dimension of velocity. It should be noted that the second 
equation in (T49|) implies that z(t) is of constant sign for t > 0. The sign of the constant cq in 
Eq. (1501) should be chosen to be the same as the sign of the initial value w, so that z(t)/co > 
for t > 0. It should also be noted that the absolute value of the dimensional constant Co is of 
no consequence for the equations of motion. The difference of two Lagrangians corresponding 
to two different constants c and c , respectively, is 



mu log — —Vi 2 + z 2 . (51) 
V Co J dt 

Since the difference (IBTj) is a time derivative, the Lagrangians corresponding to the different 
constants Co and c' are equivalent. 
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The second Lagrangian is, 



myz 2 



muzz 2 



2c cqVWT^ 2 " 
where Co is again a constant with the dimension of velocity. 



(52) 



In the Lagrangians above, the mass m occurs only as a multiplicative factor. It has been 
included for convenience in order to keep track of the dimensions in the Lagrangians in question. 

The construction of the Lagrangians (1501 and fl52|) using the methods developed in Douglas 
paper [TTJ is a task which involves lengthy calculations, which we will not present here. How- 
ever one can easily verify that the curious Lagrangians Lj and Ln above do give rise to the 
system of equations ( l4~9l) . We will consider the Lagrangian Li above only, leaving the detailed 
consideration of the second Lagrangian Ln to the interested reader. 



From Eq. ( 1501 follows readily that 



and 

whence the Euler equation 

implies that 
Likewise, 

and 



dLj_ 
dy 

<9Lj 
dy 



0. 



mz, 



dLj_ 

dy 



dLj 

dy _ 



0. 



dLj_ 

dz 



dLj 

dz 



my 



m£ 2 uz 
(I 2 + z 2 ) 



muz 



log 



1 + log 



Inserting the expressions fl57|) and fl58l) in the appropriate Euler equation 



«9L 7 

dz 



' dLj 

dz 



one obtains the expression 



'UZ 



+ 



uz 



(t 2 + z 2 )'i (i 2 + z 2 )2Z 
Using the equation fl56|) in the expression fl60|) above, one ends up with the equation 



y = t 2 uz{t 2 + z 2 y 



(53) 

(54) 

(55) 
(56) 
(57) 

(58) 

(59) 
(60) 

(61) 
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This demonstrates that the equations of motion (j49H are equivalent to the Euler equations with 
the Lagrangian Lj given in Eq. fl50|) . 

For completeness we also record the Hamiltonian Hj which corresponds to the Lagrangian 
fl50|) . The canonical momenta are defined in the standard manner, 



and 



Pz 



dz 



Py 



my 



dy 



mz, 



muz 



1 + log - 

c 



(62) 



(63) 



The equations f|62|) and (163!) can be solved for the velocities (y, z) in terms of the canonical 
momenta (p y ,p z ), 

(64) 

l + logf-^Ol. (65) 



Py 
m 



and 



Pz 

m 



uz 



This leads to the following Hamiltonian Hi, 



Hi = —PyPz + 



uz 



m 



fPy lo § 



( Py 



\mco 



(66) 



The Hamiltonian equations involving the Hamiltonian Hi in Eq. (1661) are as follows, 



dHi 

Spy 



-Pz + 



uz 



m 



z :- 



Py - 



and 



Pz 



dHi 

dz 



dHi 

dp z in 
_dHi 
dy 

£ 2 u 



:[l+l0g 

Py, 
0. 

-Py log 



Py 
mco 



Py 
mc 



(67) 

(68) 
(69) 

(70) 



(£ 2 + z 2 )i 

Combining Eqns. (IBTl) and (!7D|) one obtains the first equation in the two-dimensional system 
). Likewise, combining Eqns. (16"8"|) and (|6"9"|) one obtains the second equation in 



Using the Lagrangian Ln given above as a starting point, one obtains an alternative canon- 
ical Hamiltonian formulation of the two-dimensional system of equations (|49|) . The canonical 
momenta p y and p z are now related to the velocities y and z as follows, 



Py 



mz 



2c ' 



(71) 
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and 

\2mp, 

Pz = \ ■ 

V c o 

The Hamiltonian Hjj is 

l 2c p y 2uz 

H " = y— p ' + WTW p ~- (73) 

It is readily verified that the canonical equations with the Hamiltonian (173|) likewise reproduce 
the equations (149]) . 



2uz 



V£ 2 + z 2 



(72) 



5 Concluding remarks 

We consider the principle of stationary action for autonomous mechanical systems. The prin- 
ciple of stationary action, Hamilton's principle, can be generalized to systems with holonomic 
constraints by including the constraints in the variational procedure by means of the multiplica- 
tion rule in the calculus of variations. The corresponding procedure in the case of non-holonomic 
constraints leads to equations which are not, in general, identical to the correct equations of 
motion. This fact has been known for some time. 

It has been shown, in the case of a particular three-dimensional autonomous non-holonomic 
system, that a variational principle exists; the equations of motion for the system can be reduced 
to a two-dimensional system of autonomous equations which admit a variational formulation 
and also a canonical Hamiltonian formulation. Explicit expressions for two nonequivalent La- 
grangians together with the corresponding Hamiltonians are given. The example shows that 
one can not a priori rule out the existence of variational principles and Hamiltonians for non- 
holonomic systems, however artificial. 
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